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It was recently found that thermodynamic anomahes which arise in the Casimir effect between 
metals described by the Drude model can be attributed to the interaction of fluctuating Foucault 
(or eddy) currents [Phys. Rev. Lett. 103, 130405 (2009)]. We show explicitly that the two lead- 
ing terms of the low-temperature correction to the Casimir free energy of interaction between two 
plates, are identical to those pertaining to the Foucault current interaction alone, up to a correction 
which is very small for good metals. Moreover, a mode density along real frequencies is introduced, 
showing that the Casimir free energy, as given by the Lifshitz theory, separates in a natural manner 
in contributions from eddy currents and propagating cavity modes, respectively. The latter have 
long been known to be of little importance to the low-temperature Casimir anomalies. This convinc- 
ingly demonstrates that eddy current modes are responsible for the large temperature correction 
to the Casimir effect between Drude metals, predicted by the Lifshitz theory, but not observed in 
experiments. 

PACS numbers: ll.lO.Wx - Finite-temperature field theory, 42.50.Nn - Quantum optical phenomena in ab- 
sorbing, amplifying, dispersive and conducting media, 05. 40. -a - Fluctuation phenomena, random processes, 
noise, and Brownian motion, 42.50.Lc - Quantum fluctuations, quantum noise, and quantum jumps 



I. INTRODUCTION 



For a decade the finite-temperature correction to the 
Casimir force [l| between parallel metal plates has been a 
topic of intense investigation and debate. Describing the 
metals by a standard Drude model 
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where Q is the plasma frequency and where the relax- 
ation frequency 7 does not vanish at T = 0, the Lif- 
shitz theory implies that the temperature dependence is 
considerably different from perfect reflectors [^: a signif- 
icant thermal contribution is predicted already at dis- 
tances shorter than the Wien wavelength hc/{2TrkBT), 
on the one hand, and there is a difference of a factor ^ in 
the large-distance limit, on the other. Puzzlingly, such 
a large temperature dependence is not found in recent 
precision experiments at Purdue 13| . For reviews of the 
thermal debate around the Casimir effect, cf. [J, M, and 
references therein. 

The thermodynamics of the Casimir effect has been of 
particular interest in this context. For metals described 
by (jl.ll) , the Gibbs-Helmholtz free energy of the Casimir 
interaction is non-monotonous as a function of temper- 
ature, leading to a negative Casimir entropy in a large 
temperature range [41] . Moreover, if 7 vanishes faster than 
linear as the temperature T — ^ 0, the Casimir entropy 
remains nonzero in this limit; this was argued to violate 
Nernst's theorem, the third law of thermodynamics Q. 



Recently two of the present authors investigated the 
contribution to the Casimir force from Johnson-Nyquist 
noise, focusing on specific solutions of the Maxwell equa- 
tions for the two-plate set-up, namely purely dissipative 
(i.e., overdamped) modes which are physically Foucault 
current or 'eddy current' modes [7|. [A related investi- 
gation with a simplified model is due to Bimonte [8].] It 
was shown that the eddy current contribution alone ac- 
counts for the apparently anomalous thermodynamics of 
the Casimir effect. The non- vanishing entropy that ap- 
pears when first 7 — > and then T — >■ (taken in this 
order), is due to an infinite degeneracy of quasi-static 
Foucault current states, a glass-like situation for which 
the Nernst theorem does not apply Q. The situation is 
closely analogous to that of a free particle coupled to a 
heat bath [lOj, which is essentially in its high tempera- 
ture limit for any nonzero temperature when no damping 
is present, and for which the Nernst theorem is satisfied 
for a fixed friction rate [lli | . The apparent thermodynam- 
ical anomaly in the Casimir context was investigated in 
detail in Refs. |12h14| . It is now established that the 
Lifshitz theory gives a low-temperature expansion of the 
Casimir free energy between two Drude metals in the 
form M, EM^ 



(1.2) 



where the free energies are split into 

F{T) = Fo + AFiT), 



(1.3) 



Fq being the zero temperature value. The derivation of 
these results in Refs. |16l - [lg |. starting from a Matsubara 
sum, is quite tricky, see Sec. 3 of Ref. [i3|. They were 
confirmed independently using a scattering approach by 
Ingold and collaborators |10| . 

In this paper we go one step further in showing how the 
behavior of the Casimir effect between good Drude met- 
als is dictated entirely by the contribution from Foucault 
current modes. We introduce a free energy of interaction 
between Foucault current modes in two Drude plates, 
and find it to have the same form at low temperatures: 



AFd{T) = f^^h' + fl^^'^T'/' 



(1.4) 



We use throughout the subscript D to denote the eddy 
current (or diffusive modes) contribution to the Casimir 

(2) 

interaction. We are able to calculate the coefficients fj^ 
and /]j and, quite remarkably, find 
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The calculations are based on an analysis of the zeros and 
branch cuts of the dispersion function for the Casimir 
energy, similar to previous work based on the argument 
principle [19l |20[ . This analysis permits us to identify a 
density of states (DOS) for both the Foucault-current in- 
teraction and the full electromagnetic Casimir interaction 
within the Lifshitz theory. This method reveals a close 
relationship between the two interactions, and yields the 
results (jl.Sp in a fairly simple way, including the correc- 
tion term to /'^•' of order 7^/ri^ which we calculate in 
the limit of good conductors. 

The low-temperature expansion is valid on a tempera- 
ture scale lower than 
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where D is the diffusion coefficient of Foucault currents 
[21| and L the distance between the plates. This scale 
('Thouless energy' [22J) has been identified in previous 
work [i^ [2J, [2J] and emerges naturally when spatially 
diffusive modes in two half-spaces are coupled by elec- 
tromagnetic fields across a gap of width i [3]- It cor- 
responds to a temperature around 20 K for L — 100 nm 
and the conductivity of gold at room temperature. We 
shall refer frequently to this parameter in the following. 
The paper is structured as follows: in section |lll we 
introduce a general scheme for calculating the low tem- 
perature expansion of the Gibbs-Helmholtz free energy 
from DOS functions, and recapitulate the DOS for the 
Casimir-Lifshitz and Foucault current interactions, re- 
spectively. We use a method of contour integration to 
derive a relation between the DOS of the two types of 
interaction. This provides an intuitive tool for calculat- 
ing the desired expansion coefficients f^"^^ and /(^/^) in 
section Hill first to leading order in the small parameter 



7/ri, then the correction term. Various mathematical 
results are collected in the appendixes. 

Throughout the calculation we assume the material be 
described by (11.11) , and let ft = fee = 1 . We shall use the 
terms eddy current and Foucault current interchangably. 



II. MODE DENSITIES 

A. Introduction 

The Gibbs-Helmholtz free energy F for a system with 
a continuous distribution of bosonic normal modes is re- 
lated to the DOS p{uj) (modes per angular frequency) by 
the relation 

00 

AF{T) = T fdu p{lo) log(l - e--^/^) 
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where A^(a;) is the integrated mode density: 
p(a;) = -d^M{ijj). 



(2.1) 



(2.2) 



(We fix the integration constant with A^(0) = 0.) The 
mode density p{u)) (per angular frequency) specifies the 
physical system. (Note the difference to the density of 
states per unit energy introduced in Ref.[25|.) In the low- 
temperature limit, the exponential confines the integrand 
to small values of w, and we can expand A^(aj) in pow- 
ers of w [see also Ref. [2g|]. Integrating termwise, each 
power Lo'^ of the expansion yields a contribution ^ j'^'+i 
according to 



dLOUj" 
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(2.3) 



This method is the real-frequency analog of the method 
laid out in [lj| and used in [17| where Matsubara sums 
were expanded at low temperatures. The exponen- 
tial cutoff from the temperature dependence makes the 
procedure considerably more straightforward here, since 
standard methods of asymptotic expansion are applica- 
ble. 

For the two-plate geometry, i^ is a free energy per area 
and also depends on their separation L, with the cor- 
responding pressure given by p = —OF/dL. The low- 
frequency expansion of the mode density Md^uj) for dif- 
fusive modes is found to be of the form 



A^d(w) 
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where the inverse diffusion constant l/D conveniently 
provides the physical units, and the Thouless frequency 
D/L^ gives the relevant frequency scale. The coefficients 
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ra)j' ' are dimensionless, the first of which relates 
quite obviously to the static value of the mode density 
[see dU)]: 



.(1) 



Pd(0) = - 



D 



(2.5) 



Applying the identities ()2.3|) . we get the desired free en- 
ergy expansion 



.(1) 



AFd{T) = a2)^T^ 



+ 



0FC(i)^(3/2)T^^^(^3) (2.6) 
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where C(2) = tt^/G and C(|) ~ 1.34149. As we calculate 
in Section Hill below 

(1) _ 21og2-l ^ xj^/n)^ 
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where an expansion for good conductors (7 <C fi) has 
been performed, with corrections to m]^ appearing at 
the order 0^{'-f/fl). The plasma penetration depth is de- 
fined as A = c/H.. Note that the limit L — ?> 00 cannot 
be applied here, since it conflicts with the small param- 
eter TL^/D in the expansion [Eg. dl.Gp ]: this is why the 
scaling with L in the third term on the right hand side 
of of Eq. (|2.6l) is not unphysical. In the limit L — > 0, 
AFd is nonzero and finite: this can be attributed to the 
change in the bulk self-energy of the electromagnetic ex- 
citations of the metallic medium, as a pair of surfaces 
is introduced (the 'cleavage energy' discussed by Barton 
[27|). The surfaces introduce boundary conditions for the 
fluctuating electromagnetic modes (eddy currents in this 
case), leading to a change in energy per area with respect 
to a uniform bulk medium. 

We identify in the two following Sections the mode den- 
sities pchi^) and pd{^) that determine, respectively, the 
free energy due to all modes and due to diffusive modes. 
The former quantity is calculated within the Lifshitz the- 
ory for the Casimir effect [2a| . 



B. All modes: Lifshitz mode density 

Let us recall that the mode density pcl(w) counts how 
the mode number at a given frequency lo for two half- 
spaces at separation L differs from the situation of two 
plates at inflnite distance. The Lifshitz formula for the 
Casimir free energy [29| can be written in the form of 
Eq. (|2.1[) so that the following form of the mode density 
can be read off 



Here and henceforth, let z denote a complex frequency. 
The "dispersion function" 'D{z) is given by the integral 

V{z) ^Y.^" 

cr— p,s 
cr=p,s p 



Here, k = -^//c^ — z'^/c^, cr is a polarization index, and 
L the cavity width. In the following, we only consider 
the s- (or TE-) polarization and drop the polarization 
label. The reflection coefficient becomes (using the Drude 
dielectric function ()1.1|) ') 



r{K.,z) = rs(K, z) = 



K- J K^ + l^liz) 



K + , /k2 + Kjiz) 
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K^[Z) = — 
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(2.10b) 



All square roots are chosen here with positive real 
part; this implies in particular that ImK < and 
Im^/K^ ^ K'i{z) < for z in the upper half-plane. 



C. Eddy current modes 

The dispersion function Viz) is analytic in the upper 
half-plane. When it is analytically continued, singulari- 
ties appear on the real axis and in the lower half-plane: 
branch points where the argument of the logarithm in 
Ea. (|2.9p vanishes, and branch cuts from the square roots 
involved in the reflection coefficients (see FigH]). These 
singularities are related to the electromagnetic resonance 
frequencies of the two-plate setup that determine the Lif- 
shitz free energy from the argument principle [9|, Il2l |20| . 
They also provide a physically well-motivated way to iso- 
late the contribution of a particular class of modes to the 
Casimir interaction. 

The eddy current (diffusive) modes, for fixed k, are 
identified as a branch cut of 'D{z) along the negative 
imaginary frequency axis (see figure [ij, z ~ —i^o ■ ■ ■ ^ ij 
(^0 — £,o{k) is defined below). This branch cut is an 
example of a dispersion function that is not real on the 
imaginary frequency axis, in distinction to the familiar 
behavior in the upper half-plane. Indeed, one can con- 
firm from the macroscopic Maxwell equations that purely 
imaginary eigenfrequencies appear in a planar cavity of 
two half-spaces described by the Drude dielectric func- 
tion 121 . As is well-known in scattering theory (see, e.g. 
Ref.[20,l31|), the branch cut can be interpreted as a dense 
coalescence of discrete modes, and the relevant quantity 
is a mode density given by 



Pcl{^) 



-lindu,V{z^u + iO). 



(2.8) 



Aiz3(e) = -Ima^P(z = -ie + 0), 0<e<7- (2.11) 



Imz 



-iy 



OM-iO 



(a) 



z = ck 
f ... 



.>-Rez 



O)B(k) 



Imz 



(b) 



Imz 



(C) 



■vRez 



' : ] (eddy 



(all modes) Rez 



lTJ current C+ (propagating modes) 

modes) 



FIG. 1: (Color online) (a) Complex eigenfrequencies in the parallel plate geometry, for a fixed wavevector k (not to scale). 
The thick dotted lines show branch cuts (three-dimensional mode continuum), crosses mark poles (discrete frequencies). The 
circular dot is a pole of the first factor in Ea. (l2.18|l . The pole structure is symmetric with respect to the imaginary axis, 
according to Ea. (|2.f9|) . The frequency Ldsik) marks the transition from discrete cavity modes to a continuum of bulk modes 
(sec 30] for details), (b) Integration path: Cd around the eddy current continuum, Cui around the pole at 2: = ij + iO, C+ 
around cavity and propagating modes. The corresponding paths in the left half-plane are denoted C-ui and C_, respectively, 
(c) Integration path that encircles all modes, as relevant for the Lifshitz theory. Closing this contour in the upper half-plane, 
one gets the residue from the pole 2: = a; -f iO in the upper half-plane. 



The dispersion function is evaluated here to the right of 
the branch cut. Continuing analytically from the upper 
half-plane, we find that k is mainly real, while K-y(— i^ + 
0) = —ikj(S.) becomes mainly imaginary with 



k^iO 
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As a consequence, [k^ — fc^(C)]^^^ moves to the (negative) 
imaginary axis if k is small enough; more precisely, we 
require 



< K < fc^(^) 
This is equivalent to 

Uk) < e < 7 

where the lower bound S.oik) solves 

(7-eo)(c'A:2 + eg) = r)%. 



(2.13) 



(2.14) 



(2.15) 



We note the limiting behavior £,Q{k) w Dk^ as fc ^- 
where D is the diffusion coefficient of Ea. (jl.6l) . In the 
range (|2.14l) . the reflection coefficient becomes the uni- 
tary number 



r(K,-i^ + 0) 



K + iy'fc2(^)-K2 

>,-UkUO-^' 



(2.16) 



where the sign of the square root applies on the right 
side of the branch cut and follows by carefully evaluat- 
ing the imaginary parts of k and kj{^). For imaginary 
frequencies outside the range (|2.14p . the reflection coef- 
ficient is real (— 1 < Tg < 0), and the eddy current mode 
density (|2.1ip vanishes. 

After integrating over k, one gets a mode density fj.^ (^) 
that is nonzero in the range < ^ < 7. Finally, the 



density for eddy current modes poi^) at real frequencies 
is defined by associating to each overdamped mode z = 
— i^ a Lorentzian spectrum centered at zero frequency 
whose width is ~ Imz. Referring to Ref.iTi] for details, 
we get 



Pn(w) 



dC e 



T^J■D{0■ 



(2.17) 



D. Contour integral representation 

In this section we derive a contour integral represention 
for the mode densities of the full Casimir-Lifshitz interac- 
tion and of the eddy current contribution. This demon- 
strates a simple relation between pchi^^) and pd[ui). We 
thus prepare the low-frequency analysis we perform in 
Sec lIIII focusing on the particular case of a good Drude 
conductor (i.e., 7 ^ fi). 

It is easy to see from the sign flip of the root involv- 
ing K^{C) in Ea. (|2.10p . that the dispersion function I?(z) 
jumps and changes into its complex conjugate across the 
branch cut z = ... — 17. This jump defines the eddy 
current DOS in Eq. (|2.1ip . The latter can thus be written 
as a contour integral in the complex plane. 



PD[i^) 



dz 
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Cd 



where the path Cd encircles the cut on the negative 
imaginary axis in the positive sense as shown in Fig[TJb). 
Now, shifting the contour towards infinity, we encounter 
the poles at z = iw from Ea. (|2.18p and other singulari- 
ties (poles and branch cuts) of dz'D{z). The behavior of 
the exponential e^^^L Jq-,. j^;] _-, 00 makes I'(z) vanish 
at infinity. Hence we conclude that the integral around 
Cd is equal to the negative residues of the poles z = ±uj 



minus integrals over the contours C± in Fig[TJb) that 
encircle the singularities near the left and right half of 
the real axis. We use here the link between the disper- 
sion function and the response (or Green) function of the 
two-plate cavity [13, 1201 ^^^^ entails the symmetry rela- 
tion 



V{-z*) ^V*{z) 



(2.19) 



As a consequence, complex mode frequencies and branch 
cuts appear in pairs on opposite sides of the imaginary 
axis. (In FigHJ only the right half is shown.) 

The residues at the poles are easily calculated from the 
contours C±u in FigdKb): 



£l^^^^^(") ^ ^ (^^^^^^ + 5.2?(-..)) 



C^+C- 



-lmd^V{uj) = /9cl(w) 



(2.20) 



We thus recover the mode density for the Lifshitz theory 
as one term in the eddy current DOS. This is actually 
not surprising, since pcl(w) can be written as a simi- 
lar contour integral as Eq. (|2.18p . but evaluated along a 
contour just above the real axis [FiglTJc)] and closed at 
infinity in the lower half plane. This contour encircles all 
singularities of the dispersion function as it should, since 
the Lifshitz theory accounts for all modes. If this contour 
is shifted through infinity into the upper half-plane, only 
the two residues calculated in Eg. (12.201) contribute since 
the dispersion function is analytic inside the contour. 

In conclusion, we can write the following splitting of 
the mode density for the Casimir effect 



Pcl(w) = /Od(w) +P±(w) 



(2.21) 



where the last term gives the contribution of modes near 
the real axis (contour C+ in Figlljb), and corresponding 
C_ in the left half-plane). By continuity with the lim- 
iting case 7 — >■ 0, we can identify the latter modes with 
propagating modes in the vacuum cavity, in the bulk, or 
with electromagnetic surface modes (for example, surface 
plasmons that appear in the TM-polarization). We shall 
see in the next section that for nonzero, but small 7 ^ J7, 
the mode density p±(w) becomes small at low frequencies 
(w < D/L^,7, c/L), so that in this range, the full elec- 
tromagnetic DOS pcl(w) nearly coincides with the eddy 
current DOS poit^)- 



III. LOW-FREQUENCY EXPANSION 

We calculate now the small uj expansion of the den- 
sity of states for eddy current modes. According to Eq. 
p.21[) . we start with the full Casimir-Lifshitz interac- 
tion and discuss then the differences between the two. 
We begin with a general estimate of the scaling for good 
conductors. 



A. Scaling for weak damping 



The analysis in the complex plane, as illustrated in 
FiglH suggests that the density of diffusive modes is 
concentrated in a range ~ 7 near zero frequency. An- 
ticipating from the analysis below a total number of 
~ 1/L^ modes per unit area, one gets for 7 — >■ a scaling 
P£)(w) « 1/{jL'^)9{u}/^,L/X) where the frequency ap- 
pears in the function 9 only in the dimensionless form uj/j 
(and similarly for the distance L/X = ilL/c). A different 
behavior emerges for propagating modes (inside the con- 
tour C+ in FigH]): they move onto the real axis for small 
7 and contribute to p±{uj) in the range w ^ c/L, fl. Their 
contribution at much lower frequencies that interests us 
here, is proportional to their imaginary part and there- 
fore scales linearly with 7. This observation shall provide 
us with a simple rule to identify the respective contribu- 
tions of eddy current and propagating modes in the full 
(Lifshitz) mode density, Sec lIII Bl Note that we consider 
here the case of a fixed (temperature-independent) scat- 
tering rate 7. 

Some further corroboration of these rough estimates 
is desirable. Let us consider for the simplicity of argu- 
ment that the dispersion function dz'D{z) has only dis- 
crete poles ujn{k) in the lower half-plane, labelled by the 
momentum quantum number k. This can be achieved by 
enclosing the system in a finite box [12, 20]. One recov- 
ers the branch cuts by taking the box size to infinity [SlJ . 
From the symmetry relation p.l9p , the poles occur either 
on the imaginary axis (as for diffusive modes) or pairwise 
in the lower left and right quadrants (as for propagating 
modes). The two terms pd{'^) and p±(w) collect these 
poles, respectively. 

We make the replacement '^^Pd{C) '^ 

<^^k/{2TTfY.n&eddy and find that the DOS for dif- 
fusive modes Pd{^) can be written in the following 
scaling form 
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kdk 
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{Uk)h? 



where the limit of the mode branches for two separate 
plates (L — )■ cx)) is subtracted. We have already seen 
that the mode frequencies satisfy S,n{k) < 7. As a conse- 
quence, the integral tends toward a finite limit as 7 — >■ 0, 
Pd{'^) depends only on the scaled frequency aj/7 and is 
proportional to the scaling factor I/7. This implies in 
particular that the integral over the diffusive mode den- 
sity can give a nonzero contribution even as 7 — >■ 0. We 
confirm these results in Eq. (|3.9l) below. 

The density of propagating modes p±(w) shows a dif- 
ferent scaling with 7. With the same re-writing, the 
contours C+ and C_ collect the modes away from the 



imaginary axis and lead to the representation 

dz z 



P±{^) 



C++C- 27r z2 - a;2 



°°kdk 



27r2 



Im 



UJn{k) 



Y 

^^-^ ^2 - a;2 (fc) 

n(Eprop 



(3.1) 



L—^oo 



where in the second hne, we represent the modes by the 
poles in the lower right quadrant. Now, the imaginary 
part of the eigenfrequency Wn(fc) is negative and, for a 
small Drude scattering rate, of the order 7. A typical 
scale for its real part, on the other hand, is the low- 
est cavity eigenfrequency c/ L or the plasma frequency 
fi. Although we do not need them here, recall that the 
surface plasmon modes appear at « 17/ v2 — i7/2 for 
k ^ 0/c [33]. For an estimate of p±(a;), we focus on fre- 
quencies uj much smaller than the real part, lo <^ c/L, f2 
and take 1/L to estimate the relevant wave vectors. This 
gives 







p±H^o( 



f72L2' 
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(3.2) 



as we confirm in Eq. (|3.17p below. This small "tail" of the 
mode density near zero frequency can be understood from 
the broadening of the discrete modes due to damping (a 
(5-peak becomes similar to a Lorentzian, see also Ref . |25| ) . 
In particular, it vanishes in the limit 7—^0 where p±{'jj) 
goes over into the mode density of the plasma model 
which scales proportional to ijp . 

To summarize this estimate, we expect from Eq. p.2ip 
that as 7 — >■ 0, the low-frequency mode density for the 
Casimir-Foucault interaction and for the full Casimir in- 
teraction coincide in order I/7, with a small difference 
~ 7 arising from propagating modes. These contribu- 
tions are calculated in the following sections. We are 
thus able to check our approach against the free energy 
expansion of Refs.fla, Il7 |. 



B. Lowest order: Lifshitz theory 

Let us calculate the coefficients m}^ and nij^ defined 
in Eq. (|2.4p . starting with the leading order in the small 
parameter 7/12 which, as we have just seen, is provided 
by the full Lifshitz theory. It is convenient to work with 
the integrated mode density which from Eq. 
read off as 



Mci.{<^) =Im2?(w + iO). 



(3.3) 



We therefore start by analyzing in detail the behavior of 
^{z) for small frequencies \z\ ^ 7. 

The fc- integral in Ea. (l2.9p is re- written in terms of a 
real variable 2/ > defined by k = yK^{z). This is equiv- 
alent to a shift of the integration path in the complex 
K-plane along a more convenient direction: one still has 
convergence from the exponential exp(— 2kL) because 
Re Kry{z) > (keeping clear of the branch cut for z on the 



negative imaginary axis) . The reflection coefficient (|2.10p 
becomes real along this direction and independent of z: 



r{y) 



y - \Jy^ + 1 _ 1 

V + VFTT {y + v/y" + 1)2 



We get 
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(3.4) 



(3.5) 



where the lower limit is given by the complex number 



X{z) 



CK-y(z) 



\z z + Yy 



(3.6) 



For \z\ < 7, we have \x{z)\ < V2{-f/fl) <C 1 for a good 
conductor, and to leading order, we can replace the lower 
limit in Eq. p.Sp by zero. This defines l^oiz) and via 
Ea. p.3p . Mcl,q{^)- We write McLni'^) for the error 
(i.e., the integral from to x(z)) and calculate it in 
Eq. ipm)) below. 

By inspection, 'Dq{z) depends on z only via the func- 
tion K^{z) that can be written as 



K^{z) 



i^iizh) = -^ 



z/7 



zh- 



(3.7) 



involving the scaled quantity 2/7. From the reflec- 
tion coefficient r[y), the relevant integration domain is 
< y ^ 1. We can therefore expand the exponential in 
Ea. (l3.5p provided Ki{z/"f)L <C 1. This yields the condi- 
tion (|z^|/Cl)^^^ ^ 1 where ^^ is the Thouless frequency 
introduced in Eq. (jl.6p . Expanding to the first order in 
this small parameter, we get [D is the diffusion coeffi- 
cient) 



-Doiz) 



\D 



/g|log[l-.^(,)] 



(3.8) 






z \i/2 fy'^Ay 2r2(y) 



a^ J 2^2 l_r2(y) 





O^z/^, 



where powers z and z'^" have appeared. The integrals 
can be solved exactly (see Appendix Ej) . and we get from 
Eq. p.3p the following approximation to the Lifshitz in- 
tegrated DOS 



Mci.fl{uj) 



21og2- 1 cj 



87r2 



D 



LV2 /U}\3/2 

247r2 Vz?/ 



(3.9) 



valid for cj <C S,l,1- This proves the first term in 
Eas. (|2.7p . It is clear from this calculation (a power se- 
ries in (w/^l)^^^) that these results cannot be applied for 
7 — >■ at fixed a; > 0. In other words, the limits 7 — >■ 
and w — >■ do not commute. For a discussion, see Refs. 

El [a. 



Calculate now the small correction A^cl,7('^) from the 
lower integration limit in Ea. (l3.5p . This arises between 
the boundaries y — and y — xi^)- Recall that in the 
limit 7 ^ ri, we have |x('2)| ^ 1 and expand the inte- 
grand for 2/ ^ 1. This gives 

Xiz) 

McLni^) - -lmK'^{z)J gl log [2y(2 + k,{z)L)] 

(3.10) 



J_uP 



one half the mode density for the so-called plasma model 
where 7 = is taken from the outset. Notably, this term 
gives a contribution to the free energy proportional to 
T^, which exactly coincides with the expression given in 
[llli^l. Note that the term scaling with Eg. (1X2)) has not 
appeared in the full (Lifshitz) mode density. We outline 
an interpretation in Sec lIVI 



C. Eddy current modes 

We now address the density of eddy current modes 
alone that involves according to Eq. (|2.17p an integral 
along the branch cut of 'D(z) on the imaginary axis. It is 
again convenient to work out the integrated mode den- 
sity Md{^)- a partial integration leads to the integral 
representation 



7Wz3(w) = - 



dC 



TT ^2 



jMd{0- 



(3.11) 



where Mo(0 is the integrated mode density along the 
branch cut. By changing the momentum variable from k 
to K, this function can be written as 



k-rii) 



Md{0 



ndK 



2tt 



2lmlog[l 



r'U, 



-iC)e- 



-2kLi 



(3.12) 



«/c 



where the integrand is zero above the upper integration 
limit kj{(^) that was defined in (|2.12l) . 

The limiting behavior of this expression for a good con- 
ductor can be worked out similar to Eq. p.Sp . Writing 
the integral in terms of x = ^/7, we see that A^^l"^) 
[Ea. (|3.1ip ] depends on the scaled frequency uj/-f. The 
upper integration limit takes a form similar to Eq. ()3.7|) . 



fc^(^) = ki{x) 



1-x' 



(3.13) 



while the lower one, ^/c — xj/c, can be taken as small 
compared to the typical values k ^ 1/L and k ^ ki{x) 
that appear in the integrand. 

This motivates again a splitting of Af£)(^) in two terms, 
a first one where the lower boundary in Eq. (|3.12p is taken 



as zero, and a correction, similar to what we did after 
Eg. p.Sp . The two terms produce a split of the mode 
density p.ll|) into 



where the first term can be written as 



MdA^) 



dx w/7 

n x'^ + {oj/"iY 


ki_{x) 



(3.14a) 



(3.14b) 



j ^Inilog[l-r2(^,-ix7)e-2«^] 



Here, we have succeeded in removing from the integrand 
all dependence on 7, except for the frequency scaling. 
The second term, MD,-y{^)i is discussed in Sec lIIIDl 
Eq. p.l7p . This term is related to the correction pro- 
portional to 7 identified in Sec lIII Al the only difference 
being that we are dealing here with integrated mode den- 
sities. The expression A^£)^o(ci;) [Ea. (|3.14bl) ] is nonzero in 
the limit 7 — >■ 0, except for the appearance of the scaled 
frequency w/7. Therefore, this term corresponds to the 
(differential) mode density scaling with I/7 of Sec lIII Al 
Since we know from Eq. (|2.2ip that the leading orders for 
7 — > coincide for the diffusive modes and the Lifshitz 
theory, we can conclude 



A^D,o(w) = AlcL.o(w) 



(3.15) 



provided the frequency lo is below the range where other 
(propagating) modes appear that are not contained in 
A^£)(a;). The identity (|3.15p is checked by a direct cal- 
culation in Appendix [Bl 



D. Damping correction of eddy current modes 

We now show that one gets for good conductors the 
second term, of relative order (7/ri)^, in the coefficient 

nij^ [Eq. ()2.7ap ]. It arises from the correction TM^ ,y(w) 
to the diffusive mode density. It is interesting that this 
shows a scaling ^ 7W, in distinction to the correction in 
the Lifshitz theory, Eq. ((3?T0)) . 

The second term in Eq. (|3.14ap . M.r).i{'^)i is of the 
same form as Eq. (|3.14bp . with the upper limit ki{x) re- 
placed by ^x/c. For good conductors, the upper integra- 
tion limit K < jx/c is small compared to the scale ki{x) 
[Eq. (|3.13p ] that appears in the reflection coefficient. The 
argument of the exponential is small if we take 7 <C c/ L. 
Expanding both quantities for small k, we get 



MD.-riuj) 



dx uj/"f 





^x j c 



n x'^ + {uj/jy 



(3.16) 



ndn 

'2^ 



Imlog[2K(L + 2i/fci(x)). 
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The imaginary part does not depend on k, and the in- 
tegration gives a factor i(7a;/c)^. At this stage, we can 
take the low-frequency Hmit {uj ^ 7) and are left with 



MD,ii(^) 






bJ 7 



^ Ax 



— arctan | —^ 
A 



2A l-x 



D4TT^n^2X + L' 



(3.17) 



where A = c/Q is the plasma wavelength. This yields 
the correction to rrijj appearing in Eq. (|2.7a[) . We have 
checked that A4D,-y{'^) does not contain, at the next or- 
der, a fractional power w'^/^, as found for A^d^o(w). 

We suggest the following interpretation for this cor- 
rection: it is related to the mutual influence of the two 
types of modes, overdamped and propagating waves. To 
wit, as the two slabs approach each other, the differ- 
ent mode frequencies cannot shift independently because 
taken all together, they have to satisfy a sum rule quoted 
in Ref.[i3: 



h 



y^ lmuJn{k) 



.all modes 



(3.18) 



where the notation assumes that branch cut continua 
have been discretized (see Sec lIII A|) . The eddy current 
modes play a crucial role in satisfying this sum rule. In- 
deed, any modification in the imaginary part of the prop- 
agating (cavity and bulk) modes due to a change of the 
distance L (i.e. the propagating modes leave the con- 
tinuum above the plasma frequency and become discrete 
cavity modes as the distance L is increased) is simulta- 
neously balanced by a shift in the diffusive mode density 
on the imaginary axis that extends down to —17. 

Due to the sum rule p. 181) . the small correction for 
eddy currents appears also, with opposite sign, in the 
propagating modes. For this reason, the Lifshitz mode 
density does not contain this term [see Ea. p.9p ]. and 
its next-order correction Eq. (|3.10p is independent of the 
damping rate 7. 



IV. DISCUSSION AND CONCLUSIONS 

We have calculated the low-temperature behavior of 
the interaction between two parallel half-spaces across 
a gap of width L due to low-frequency Johnson noise in 
the bulk of the conducting medium, in particular eddy or 
Foucault currents that are coupled to TE-polarized elec- 
tromagnetic fields. The interaction is calculated in orders 
T^ and T^/^ and is compared to the Casimir free energy 
within the Lifshitz theory for Drude metals. A striking 
result is uncovered: the low-temperature correction to 
the Casimir effect between parallel slabs of good Drude 
conductors is dictated entirely by the contribution from 
eddy currents, as demonstrated by the two leading order 
correction terms as T — > 0. This adds a further piece 



of support to the findings of Ref.fT"] where the unusual 
physics of the thermal Casimir effect between Drude con- 
ductors is attributed to the interaction between eddy cur- 
rents. 

Within our approach, we find small differences to the 
free energy that are of second order in the ratio scattering 
rate to plasma frequency, 'y/fl [Ea. (l2.7ap ]. This correc- 
tion refiects the mutual influence between the modes that 
are constrained by a sum rule, Eq. (|3.18p . Note the cu- 
rious fact that this makes AFu depend on L already at 
order T^, different from AFcl. Therefore the eddy cur- 
rent interaction makes a tiny contribution to the pressure 
(p — —dF/dL) quadratic in temperature. However this 
is exactly cancelled by a corresponding contribution from 
propagating modes and the resulting Casimir pressure is 
proportional to T^/^ to leading order, as Eq. (|3.9p shows. 

Let us finally emphasize the analysis of the singulari- 
ties of the Lifshitz dispersion function that we performed 
in the complex plane. This picture identifies in a natu- 
ral way the mode frequencies of the system, even in the 
presence of dissipation, and justifies a natural splitting of 
the free energy in contributions of specific types of modes. 
We gained in particular the insight that the mode density 
for the full Casimir-Lifshitz interaction is simply the sum 
of eddy current modes and of propagating cavity and bulk 
modes. The second contribution becomes small at low 
frequencies, weak damping and not too large distances 
because the complex mode frequencies are located suffi- 
ciently far away from the origin. This provides a deeper 
understanding why propagating modes are of little rel- 
evance to the temperature dependence of the Casimir- 
Lifshitz interaction between Drude metals. Indeed, this 
dependence was previously found to originate primarily 
in low- frequency evanescent modes [23l |24| . 
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Appendix A: Integrals for Lifshitz theory 



The integrals in Eq. (|3.8p can be evaluated with the 
substitution y = sinht. This simplifies the reflection co- 
efficient (|3.4p into r{y) = — e^^*. Hence 



/ 9 / N \ / , sinh 2i , , 

dyy\og{l-r\y)) = / di — - — log (l 



(Al) 



Expanding the logarithm, integrating term by term and 
evaluating the sum, we get 



!!^,„,(l_e-«, = -ii!5i^. (A2) 



For the second integral in Eq. (|3.8p . the same substitution 
gives 



dJ/i37:2(^ = ydt^^^ = -. (A3) 





Appendix B: Integrals for eddy currents 



We prove here Eq. (j3.15p : the low- frequency mode den- 
sities for eddy currents, Md.o{^)i and for aU modes, 
A^CL,o(w), coincide to leading order in 7. 

Consider Eq. ()3.14bl) for the eddy current mode density. 
We want to write this as a contour integral, similar to 
Eq. (|2.18p . around the eddy current branch cut Cd in 
FigH] Note first that the Im can be pulled in front of 
the K-integral and that integral be extended from ki(x) 
to 00. This is possible without changing the value of the 
integral if k is taken just below the real axis, the refiection 
coefficient (j2.16p getting real and smaller than unity in 
modulus. Hence, the logarithm is real, and this part of 
the integration range does not make any contribution to 
the imaginary part. 

The contour integral in the variable z — ~i'jx±0 finally 
takes a form similar to Section Pi Dl 



A^D,o(w) 



dz 



Cd 



27ra;2 



^Vd,o{z) (B1) 



and the reflection coefficient is given by Eg. (12. 101) . Note 
that to the right of the branch cut, fci(x) = iK-y(z). 

The variable change k — yK^{z) with y > now 
shows that the function Vufl^z) is indeed identical to 
the small-7 approximation to the Lifshitz dispersion func- 
tion, I?o(z), defined by setting the lower integration limit 
in Eq. (|3.5l) to zero. Note that this actually shifts the inte- 
gration path in the lower right quadrant of the complex 
K-plane: from I{,eKj{z) > 0, convergence at y — > cx) is 
secured. The refiection coefficient r{y) [Ea. (l3.4[) ] is ana- 
lytic and of modulus smaller than unity in this quadrant, 
hence the logarithm encounters no branch points. 



We still have to evaluate the integral (jBTj) and do this 
with the same technique as in Sec lIIDl Pulling the con- 
tour Cd to infinity, one encouters simple poles aX z — iw. 
In the present case, we can argue that the function 2?o(z) 
is analytic in the right half-plane and by the symmetry 
relation (|2.19p also in the left half-plane: this is due to 
the way the integration variable y keeps the wave vec- 
tor K, clear of the branch cuts that are located inside 
the contours C± [FigH]. Indeed, across these cuts ei- 



where Vd^q^z) is the integral 



ther K, or \ k^ + i^'^{z) are purely imaginary and jump in 

sign. This never happens along the path parametrized as 
K = yKy{z), as can be checked easily. Indeed, if z is in the 
right half-plane, k remains in the lower right quadrant, 
excluding the real and imaginary axes. 

As a consequence of I'd (-2^) being analytic in the left 
and right half-planes, the integral (jBip is given by the 
pole contributions 'Dq{±ijj) only. Referring to Ea. p.3p . 
we thus get the desired link to the approximated Lifshitz 
mode density 



Md,o{^) = ImVoiuj) = McL,o(^) (B3) 



00 

2?Ao(2)= /^log[l-r-2(«;,z)e-2KL]^ (22) which is Eq. (013. 

J 27r^ 
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